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Abstract Weak decay tree and annihilation - t-channel W- exchange amplitudes for 
the D° — > Kg-K^TY^ process are calculated using quasi two-body QCD factorization 
approach and unitarity constraints. Final state strong Kn and kit interactions in S, 
P and D waves are described through corresponding form factors including many 
resonances. Preliminary results compare well with the effective mass distributions of 
the Belle and BABAR Collaboration analyses. 
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1 Introduction 

Why should one study the weak decays D — > Kgn + n~? First, the recent mea- 
surements of the D°-D° mixing parameters for this self-conjugate reaction by the 
BABAR [1] and Belle [2] Collaborations could show the presence of new physics contri- 
bution beyond the standard model. Second, the Cabibbo-Kobayashi-Maskawa, CKM, 
angle 7 can be evaluated from the analyses of the — > D^K^,D° — > Kg-ir + ir~ 
decays. Third, one can learn about the final state meson-meson interactions, the me- 
son resonances decaying into different meson-meson pairs and their interferences in 
the Dalitz plot. One can also perform a partial wave analysis of decay amplitudes. In 
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addition, constraints from quasi two-body QCD factorization [3], QCDF, approach will 
allow to test theoretical models of form factors entering in the decay amplitudes. 

Quasi two-body factorization Following a program devoted to the understanding of the 
rare three-body B decays (see e.g. Ref. [4]) the presently available D° — > Kgir + -K~ 
data are analyzed in the framework of QCDF. Neglecting the small CP violation in 
K° decays, it will be assumed that \K%) = (\K°) + \K°))/y/2. The three-meson final 
states K°n + TY~ are approximated (quasi two-body approximation) as being formed by 
a meson-meson state, [K ir + ]s,P,D or [K°ir~]g p jj or [^ + ^^]s,P,D m a S,P or D wave 
created by a qq pair and the remaining meson, n~ , 7r + or K° , respectively. Amplitudes 
are derived from the weak effective Hamiltonian, -ff e // which is a superposition of left- 
handed quark current- current operators. For instance one has the operator 0\ = 
where ji = s a 7 v (l — 7 )ea = (sc)v-A and j 2 = upjv(l - 75)^/3 = (ud)v-A> a and /3 
being color indices. Applying QCDF to the part of the amplitude proportional to the 
product of the CKM quark mixing matrix elements V* s V uc i = A\ one has, 

(K° tt-tt+I H eff \D°) ~ ^| Aijax (n+\j 2 \0) ([K° w-} s ,P, D \h\D ) 

+ a 2 (#°|j 2 |0) ([* + *-]s,p,d\3i\D°) + a 2 (0\A\D°) [([K° *-] s , P , D n + \f 2 |0> 

+ (K° [* + n-] S ,P,D ba|0>] } (1) 

where Gj? = 1.166 x 10~ 5 GeV~ 2 is the Fermi coupling constant, 0^2 are effective QCD 
Wilson coefficients and j[ = (uc)y_A> j' 2 = (sd)v~A derive from via a Fierz trans- 
formation. In Eq. (fJJ it appears that the cross product of bilinear quark currents in 0\ 
factorizes, after the introduction of the vacuum state, into the product of two matrix el- 
ements. The first ones are proportional to the 7r + , K , D° decay constants f v , fx, /d° 
since (7r + |j 2 |0) = ifnp w + , (R°\ j' 2 \0) = ifxPRO, (0\fi\D°) = -if D apr>o, P^+,P K ° 
and pjjo being the n + ,K° and D° four momenta, respectively. The second ones are 
transition matrix elements or form factors as with Mi = K ,^, j = j\ t2 or j[ 2 
one has, ([M 1 M 2 ] S ,p,p \j\D Q ) = (D Q [M 1 M 2 ]s,p,D b'|0> and ([M 1 M 2 \ S ^ D M 3 \j\0) 
— ([MiM2]s,P,D I J 1-^3 )• It can be shown from field theory and using dispersion re- 
lations [5] that these form factors can be calculated exactly if one knows the D°- 
[M\M%]s p jj or M^-[MiM 2 ]s p r) strong interactions at all energies. 

2 Decay amplitudes 

Different type of amplitudes Amplitudes with c — > sud transition are proportional to 
V* s V uc l where V cs w V uc i ~ cosOq ~ 0.975, 9q being the Cabibbo angle. There are 7 
such allowed tree amplitudes: 3 for the tt^~[K "K~]s p d> 3 for the K°[k + TT~]gp r) 
and 1 for the K u)(u) — > [7r + 7r~]p by G-parity violation) final states. Amplitudes with 
c — > dus transition, proportional to sin 2 0q « (0.225) 2 , are doubly Cabibbo suppressed. 
There are 6 of them as the W meson cannot couple to the [K°-k + ]jj final state. There 
are also 7 allowed tree annihilation or t-channel VF-exchange amplitudes corresponding 
to the cu annihilation into sd and 7 doubly Cabibbo suppressed annihilation amplitudes 
from the cu annihilation into ds. Altogether, the quasi two-body QCDF approach leads 
to 27 non-zero - 13 tree and 14 annihilation - amplitudes, 
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Transition matrix elements Several meson resonances can decay into the two me- 
son final states in S, P or D wave. For the kaon-pion subsystems, the scalar reso- 
nances, ^^(SOO)^ or and ATq(1430) =I: decay into [K^TT^g, the vector resonances, 
K*(892) ± , Ki ( 1410) ± , K*( 1680) ± into [^°7r ± ] P and the tensor resonances K^UM)* 
into [RK^rj. For the pion-pion subsystems, the scalar /o(600) or a, /o(980), /o(1400) 
decay into [7r + 7r _ ] s , the vector, p(770)°, w(782), /9(1450)°, p(1700)° into [7r + 7r~]p and 
the tensor /2(1270) into [7r + 7r~]£). This leads to a rich interference pattern in the Dalitz 
plot. These resonances are used to write the three-meson transition to the vacuum as, 



{M 1 ( Pl )[M 2 (p 2 )M 3 (p 3 )] s ^ D \ j'\0) ~ G r m 2 m, (s 23 ) (MMRf^l /|o) , (2) 

where S23 = (f>2 + p 3 ) 2 ■ The vertex functions G„m 2 u 3 (s 23 ) describe the Rq p M JZ 

resonance decays into the states [M 2 M 3 ]g pp. A similar equation holds [5] for the D° 
transitions [see Eq. |l|] to two-meson states. As an example of application of Eq. ([2]| 
one can choose Mi = R° , [M 2 M 3 ] S = {tt + tt~] s and i?g +7r = /o(980), so that, 



{R°(pi)fo(P2 + P3 )\(M) V _ A \0) = -» mg ° 2 823 P D oFf M980 \m 2 D „) + 1 term, 

P D o 

(3) 

where p D o = p\ + p 2 + p 3 - The K° to /o scalar transition form factor, ^°{m 2 D0 ), 
related to the R° fo interaction at (px + Pf ) 2 = m2 jo, is a complex number to be 
fitted. The extra "1 term" gives a null contribution when multiplied, in Eq. [3] by 
(0| (uc) v _ A \D°) [see Eq. ©]. The vertex function G /o(980) ( S 2 3 ) ^ X2 Fq + *~ (s 23 ) 
where \2 is a constant fitted to data and Fq v ( s 23) is the complex pion scalar form 
factor [4] related to ^[7r + 7r~]5| (uu)y_A\0) ■ The other vertex functions for the kaon- 
pion 7i or pion-pion ,4, systems in Eq. (J2]) are described in terms of the corresponding 
form factors. The D-wave meson-meson form factors are represented by relativistic 
Breit-Wigner formulae. In the transition form factors for Rcfp^jj we choose Rg * = 
Ko(1430), flf % = JsT*(892), Rf = ,f (980) and R P * = p(770)°. 

A selected amplitude From Eqs. ([TJ) and © one obtains for the D° — ¥ Kgn + n~ S-wave 
annihilation amplitude with [7r + 7r~]5 subsystem in the final state, 

A n 2S = A x X2 0,2 fDo{m 2 K - So) Fq fo (m 2 D o) Fq + * (s ), (4) 

where s± = {p^± + Pk ) 2 i s — (Ptt+ + Ptt-) 2 - Other amplitudes can be derived 
similarly, their expression will be given in a forthcoming paper, however some explicit 
formulae can be found in Ref. [6]. 



3 Experimental D° — > K%-k+-k data 

Isobar model and unitariy The experimental Dalitz plot analyses 1 , 2 are performed 
within the isobar model to describe the final state meson-meson interactions and many 
free parameters are used (of the order of 2 per amplitude): the BABAR Collaboration 
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model relies on 43 parameters and that of Belle on 40. Amplitudes are not unitary 
neither in the 3-body channels nor in the 2-body sub-channels. Two-body unitarity is 
incorporated in the present model: unitary form factors are used in the Kn S- and 
P-wave [7] and in the nir S-wave [1] amplitudes. The branching fraction of the sum of 
these amplitudes is larger than 80% of the total D° — > Kgn + n~ branching fraction. 



4 Preliminary results and concluding remarks 

The present model has 28 free parameters, mostly unknown transition form factors and 
a minimization procedure is used to reproduce the Kgn~ , Kgn^ and 7r + 7r~ squared 
effective mass m| = s± and ttiq = sq projections of the experimental Dalitz plot 
analyses 1,2. Preliminary results are shown in Fig. [T] for the Belle Collaboration 
analysis 2 . Similar results are obtained with the BABAR analysis Q]. 

Conclusion This theoretically constrained analysis might be useful to improve the 
determination of the D°-D° mixing parameters and of the CKM angle 7. 
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